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Abstract 

In warm inflation models there is the requirement of generating large dissipative couplings of 
the inflaton with radiation, while at the same time, not de-stabilising the flatness of the inflaton 
potential due to radiative corrections. One way to achieve this without fine tuning unrelated 
couplings is by super symmetry. In this paper we show that if the inflaton and other light fields are 
Pseudo-Nambu-Goldstone Bosons then the radiative corrections to the potential are suppressed 
and the thermal corrections are small as long as the temperature is below the symmetry breaking 
scale. In such models it is possible to fulfill the contrary requirements of an inflaton potential 
which is stable under radiative corrections and the generation of a large dissipative coupling of the 
inflaton field with other light fields. We construct a warm inflation model which gives the observed 
CMB-anisotropy amplitude and spectral index where the symmetry breaking is at the GUT scale. 



PACS numbers: 98.80.Cq, ll.10.Wx, 14.80.Va 



1. INTRODUCTION 



The cosmological theory of inflation [l[ offers an explanation of the scale invariant horizon 
sized perturbations in the CMB observed by COBE [2] and WMAP A successful 

model of inflation requires the potential of the scalar field to be flat. A natural model 
for a flat potential is through Pseudo-Nambu-Goldstone bosons in a class of model called 
"natural inflation" The flatness of the inflaton potential in Natural inflation models 

is ensured by the fact that the PNGB couplings with other fields are suppressed by the 
symmetry breaking scale / and by making / large enough one can control the quantum 
corrections to the inflaton potential. One major problem of the Natural inflation models is 
that the spectral index of the inflaton perturbations is related to the symmetry breaking 
scale as n s = 1 — Mp/(8irf 2 ) and in order to be consistent with the WMAP measurement 



n s = 0.963 ± 0.014, the symmetry breaking scale has to be close to the Planck scale [10 



In reference [U| it was pointed out that in warm inflation models |12Hl4j] where there is a 
dissipative coupling between the inflaton and the radiation bath, one can have inflation with 
PNGB and the symmetry breaking scale could be lowered to the GUT scale / ~ 10 16 GeV. 
In order to have a lower symmetry breaking scale in natural inflation models the dissipative 
coupling T of inflaton with radiation must be much larger than H. 

Inflaton couplings with radiation fields through which one can generate large V are also 
expected to spoil the flatness of the inflaton potential due to radiative corrections (lol. 16 . 
A solution for getting around the opposing constraints of obtaining a large dissipative co- 
efficient but with small corrections to the inflaton potential, without fine tuning unrelated 



couplings, was obtained by invoking supersymmetry [17|, [18j. Inflaton was coupled to heavy 
catalyst fields with masses larger than the temperature of the Universe and these fields in 
turn were coupled to light fields. The evolution of the inflaton would induce light particle 
production via the heavy catalyst fields. Since these heavy catalyst fields are expected to 
populate mostly their ground state, the quantum corrections associated with them could 
be canceled in super-symmetric models 13| however the temperature dependent dissipative 
damping terms are not canceled. 

In this paper we study warm inflation caused by Pseudo-Nambu-Goldstone bosons as the 
inflaton. We study the Ante temperature corrections to the potential. We show that the 
thermal corrections to the inflaton potential are small as long as the symmetry breaking 
scale / T. We point out the possibility that if the light fields involved in the dissipative 
warm inflation are all PNGB's which arise in successive spontaneous symmetry breakings 
of the interactions of a heavy catalyst field then it is possible to generate a large dissipative 
couplings of the inflaton to radiation without destabilizing the inflaton potential. In the 
final section we show that taking the symmetry breaking near the GUT scale one can get 
warm inflation models with PNGB which^ satisfy all the observational constraints on the 
CMB-anisotropy spectrum from WMAP 
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2. WARM INFLATION WITH PNGB 



Natural inflation models use the PNGB potential of the form which arise from an explicit 
symmetry breaking, 

V{4>) = A(T) 4 (l - cos (j^y)) • (1) 

Here we assume that the spontaneous symmetry breaking scale / and the explicit symmetry 
breaking scale A are both dependent upon temperature through loop corrections. We will 
determine the finite temperature corrections to A and / taking an example of PNGB in a 
SU(N) technicolor theory where we assume fermions condensate are formed giving rise to 
Goldstone bosons. Consider a theory with N flavors of fermions with an approximate chiral 
SUl{N) x SUr(N) symmetry. The axial part of the above symmetry is broken spontaneously, 
when the fermion condensate (ipip) attains a nonzero vev. The low energy dynamics is then 
given in terms of the matrix field U(x) = exp(2i<ft(x) / /) in terms of the (N 2 — 1) Goldstone 
bosons 4>(x). Here <p(x) = <ft a (x)T a , a = 1,N 2 — 1 and T a, s are the generators of SU(N). 
We have shown in Appendix A that the temperature dependence of A and / are as follows 
(A 4 = -^W>ofEqn. (jQ)), 

a W -^(1-35^), (2) 

and 

f( T ) = f (l " N ^p) (3) 

to the leading order in T/f. 

The dynamics of inflaton field is governed by the equation 

0+(3# + r)0 + V'(0,T) =0. (4) 

Here over-dots represent derivative w.r.t t and / denotes differentiation with respect to <p 
and T is the damping term. The total energy density and pressure of the system is given by 

P = y + V(<j ) ,T) + Ts J (5) 

P=£-V(4>,T). (6) 

Here s is the entropy density of the system that is given by the derivative of the potential with 
respect to temperature. The Friedmann equation for expansion and the energy-momentum 
conservation equation are given by 

H 2 = ^P, (7) 
p + 3H(p + p) = 0. (8) 
From above equations we get the energy conservation equation for the radiation as 

Pr + 4H Pr = T0 2 . (9) 
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During warm inflation p r remains nearly constant. Assuming a slow roll of the inflaton we 
neglect <fi in the equation (j4]) and kinetic energy term in the total energy density. During 
inflation the potential energy of the inflaton field is dominant so we also can neglect the Ts 
term in the total energy density. So we get 



V 



H 2 



The slow role parameters are defined as 



3H + r ' 
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(10) 
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(12) 



and the ratio between the thermal damping and the damping due to expansion is given by 



Q 



r 

3H' 



(13) 



Here three extra slow roll parameters appear because of dependence of damping term 
and temperature dependence of the potential and dissipation coefficient T. It was pointed 
out in 19| that these all parameters should be less than 1 + Q to realize warm inflation. 
Latter, a detailed stability analysis was made by Moss et al [lo[ and they showed that 
to have successful warm inflation the parameter /3 (b in 20[) must be less than and 
|c| < 4. This implies that thermal corrections to the inflaton potential must be suppressed. 
A detailed study of the perturbation spectra in warm inflation models has been done by 
2l| for constant as well as temperature dependent T. However, we use the results of 
to calculate the power spectrum and spectral index for curvature perturbations, 



which are given as 



Pi 
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(15) 



where A = 4(1 + Q) + (Q - l)c. 

In the following sections we calculate the dissipation coefficient in a model where the 
inflaton and the light fields are PNGB fields which interact through a heavy catalyst field. 
We will compute the amplitude and spectral index of the curvature perturbation in this 
model and compare with the WMAP observations to fix the scale of the symmetry breaking. 
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3. DISSIPATION THROUGH GOLDSTONE-BOSONS 



In our model, inflaton field <p transfers its energy to the field o via intermediate fields 
X- We assume that both the inflaton and a fields are PNGB with spontaneous symmetry 
breaking scales / and f a respectively. We work in a high temperature regime where / > 
T > m x > m a , m^. 

The interaction of (p with \ is described by the Lagrangian 

£ = j 2 (#V) X 2 (16) 
and the interaction of x with a fields is described by 

C=±(d„a)(d»a) X 2 . (17) 

J a 

Such interactions are similar to Goldstone boson baryon interaction term in chiral pertur- 
bation theory (22[. Integrating by parts, interaction rtlBjl can be re-expressed as 

C=-jt (vd^x 2 + 2<pd»( PX d ll x) (18) 

Here tp = <p + <f>\, <p representing the zero mode classical field. Hence the interaction 
term (fT9j) with the classical background field can be expressed as 

C = ~J2 {<M>X 2 + + Z^faxdpx) (19) 

In slow-roll regime the <j) term can be neglected. The dissipation coefficient T can be derived 
from the effective equation of motion that for the above interaction can be given as [l8[ 



0(t) + v'W) - ^(xx) - jsiVfaxfyx) = o, (20) 

where (...) represents ensemble averages with respect to an equilibrium (quantum or thermal) 
state. These field averages can be calculated using linear response theory and (xx) can be 
written to the first order as 

(XX) = (XX)0 ~ * j d 4 x'9(t - [<j>(x')<f>(x') - j>{x)<j>(xj\ ([x(x') X (x'), x(x) X (x)}). 

(21) 

Here (xx)o represents the correlation function evaluated at initial time t. In slow- roll regime 
we can take <fi(t) as nearly constant so the effective equation of motion (|20|) becomes 



(A \)n / I ) / d 4 x'9(t - f ) (0(x') - <f>(x)) (Ix(x')x(x'), x(x) X (x)}) = 0. 



P J 



(22) 
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The non-local term (<p(x') — <f)(x)) in the integrand of the above equation can be simplified 
by Taylor expanding and the effective equation of motion becomes 



4>(t)+V(<j>) 



P 



(xx)o + i 



d*x'e(t - o (f - 1) (W)x(x'), X (x)x(x)}) 



(23) 

We see that the contribution to the dissipation coefficient T comes from the two coeffi- 
cients of the <p in Eq. ff23l . The first term is given by 



r i = —^(xx)o, 



(24) 



while the second term will be 

/ 4 



d*x'9(t - f) (f - t) ([x(x')x(x'),x(x)x(x)}} 



(25) 



In terms of the spectral functions of the x field, the first term becomes [23 



r, = --pi 



d 4 p 



iPoP x (Po,P)n x {p ) 



(26) 



where n x is the distribution function for the % fields and p x is the spectral function defined 

as as 



P x (Po,P) 



2JmS, 



(27) 



where u p = ^p 2 + rnE. Using the relation lmYj x = 2u p T x and substituting ( 1271) in ( )26l) we 

get 

d 3 p dpo 



r 1 = 4 



f 2 J (2vr) 3 2vr 



Po 



4u p r x 



n x (Po 



(28) 



For a given temperature, the po integral is dominated by the point u p , which lies close to 
the poles of the spectral function. So to evaluate p integral the integrand can be expanded 
about po = Uk and we obtain 

ri = 4 / t^%K) (29) 



that gives 



PJ (2tt)3 



T 3 
i p 



(30) 



This term is not sufficient for the strong dissipation as j « 1. To get large dissipation we 
will consider the second term ( 12 5 p which can be represented diagrammatically as Fig. [TJ As 
this is the main contribution to the dissipation coefficient so we will leave the subscript 2 in 
this term and denote it by T. 
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FIG. 1: The self energy diagram of ip which gives rise to the dissipation width T. 



Again we can write T 2 in terms of spectral functions and we get 231 ]. 

T 2 = T = YJ^j -0jlPlPx(Po,P) 2 x n x (l + n x ) (31) 

Again we can follow the similar steps as earlier to calculate T. Here after evaluating the 
energy integral we get 

T = W 4 I W?^~ x x nx{Up) (1 + nx{Up)) ' (32) 

The relaxation time t x for the \ can be calculated using the interaction ( TTT1) . Since the a 
field is also a PNGB it has derivative coupling to the intermediate field x- We have used the 
method given in 24| to calculate the relaxation time. Here since the a field is in thermal 
bath, the contribution to the relaxation time of the x ne ld will come from the thermal part 
of the free propagator of the internal lines. As it is evident from the finite temperature 
spectral function of the a field (see Eq. (A10) of 24|), the derivative term in the interaction 
can be replaced by the mass of a field using on-shell condition. Finally the relaxation time 
will be 

1 3 m A T 2 

* 1 = r < = TO (33) 



Putting this expression for T x in f[32|) and redefining the integration variable j3u p = x the 
dissipation coefficient will be 



l2 fA 



r= ^?4 T £/V^-M^^. (34) 

Now in the regime T > m x the thermal integral can be evaluated numerically and it gives 
a factor of 7. Hence the dissipation coefficient will be 

We assume that and a are generated from successive symmetry breaking such that / > 
f a . The PNGB masses m and m a are related to the symmetry breaking scales / and f a 
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respectively and some explicit symmetry breaking term in the Lagrangian like the fermions 
mass terms in a SU(N) chiral symmetry breaking model. We will assume for simplicity that 
the explicit symmetry breaking scale for both the PNGB are of the same magnitude A. Then 
from the Gell-Mann-Okubo-Rienner relation we have m? f 2 = m 2 f 2 = A 4 . Using the GOR 
relation and denoting = (fa-/ f) the expression for F in terms of can be written as 

F = ——tt—T. (36) 
3vr A 8 v ; 

There will be thermal corrections to the effective potential of </> due to the x loop that 
can be given as [25| 

m 2 6 2 / T \ , 
V* = — (37) 



8vr/ 4 \m x/ 

With the choice of the parameters required to satisfy CMB observations (see the next sec- 
tions) these corrections will be of the order of ™J/a (jl^j ~ 10~ 8 GeV 2 which are much 

smaller than the mass term of the inflaton i.e ~ 10 16 GeV 2 . 

In the next section we will use this expression in the formula (TT4"]) and ( TT5|) for the power 
spectrum amplitude and spectral index and fix the value of / which gives the correct CMB 
spectrum. 

4. OBSERVATIONAL CONSTRAINTS 

The PNGB potential parameters and warm inflation parameters can be constrained by 
using observational bounds on amplitude and spectral index of primordial perturbations. 
Using expression fl36|) for T, the inflaton equation of motion (j4]) can be rewritten as 

4>+\3H+ ^jf-T^j + V'(4>) = 0. (38) 
To determine we assume slow-roll approximation and we get 

1/3 

(39) 



3irV'A 8 



56c«/ 4 T 

We use this value of <fi in the rest of the numerical calculations and we get T ~ 10 12 GeV, 
which satisfies the conditions T 2 » V"((f>) ~ 10 17 GeV 2 and r >> H ~ 10 6 GeV. This 
justifies our assumptions used in Eq. fl39|) . 

Assuming the amount of radiation produced by dissipation is nearly equal to the radiation 
diluted due to expansion, the radiation density is given as 

AH Pr = T0 2 . (40) 



Using p r = T^g+T 4 and (135]) . fl39l the temperature of the thermal bath is given as 

(41) 



30 ; 

>4\ V13 



15 y /13 / 3-7r\ 1/13 / a 8 (y'(<f>)y 
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We have solved equation ( 1551) and H = —AttG yfi 2 + |p r j numerically using e-foldings 

iV = din a as an independent variable. To solve these equation we take the values of the 
parameters that satisfy the observational constraints on A s and n s . We take = 0.145, 
/ = 1.29 x 10 16 GeV and A = 2.45 x 10 12 GeV. The evolution of infiaton field 0, radiation 
density p r and potential V(<f>) is shown in Fig. |2]and|5J As depicted in Fig. [5J inflation ends 
when radiation density becomes equal to the potential energy of the infiaton field and we 
enter in a radiation dominated phase. With this choice of parameters, we can have sufficient 
number of e-foldings required to solve the horizon problem. 

We have three parameters c^, f and A in this model that can be constrained from obser- 
vations. In FigJU we have plotted the power spectrum and spectral index by varying / and 
A and using (155]) for T and (1411) for T. It is clear from the figure that the symmetry breaking 
scale is at the GUT scale to satisfy the observation constraints. From this figure we can 
see that for A s = 2.38 x 10" 9 and n s = 0.959, which are within la, / = 1.29 x 10 16 GeV 
and A = 2.45 x 10 12 GeV. So for <fi = 2.8 /, we get the values for other parameters using 
Eq 056]), gl]), (|32}, (HD and we find that F = 1.66 x 10 12 GeV, T = 2.09 x 10 n GeV and 
H — 1.99 x 10 6 GeV. The ratio of the dissipation coefficient and Hubble constant Q ~ 10 5 
is quite large while the thermal corrections to the infiaton potential are suppressed as T/f. 



5. CONCLUSIONS 

We have shown that it is possible to achieve natural inflation at the GUT scale / ~ 
10 16 GeV in warm inflation models if the dissipative coupling of the inflatons is large T ~ 
10 12 GeV. We compute the thermal corrections to the PNGB potential and show that the 
corrections are small as long as T C /. We also show that a large dissipative coupling 
can be achieved if the both the infiaton and a light radiation field are PNGB's which arise 
in successive spontaneous symmetry breakings and which couple to a heavy catalyst field. 
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FIG. 4: The allowed range of /(GeV) and A(GeV) using the WMAP 7 constraints on amplitude 
A-ft 2 and spectral index n s of curvature perturbations. 



Thus PNGB model of warm inflation is another class of models besides supersymmetry 
where one can naturally fulfill requirements of generating a large dissipative coupling without 
destabilizing the inflaton potential by thermal and quantum corrections. 



It has been shown recently |26|, |27J that in a model of inflation with PNGB there is a 
non-gaussianity due to the coupling of the pseudo-scalar with gauge bosons. In models such 
as ours where f « M P the non-gaussianity expected to be large and this could be observed 



in the forthcoming Planck [28[ observations. 
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Appendix: Thermal corrections to the PNGB potential 

Consider a theory with N flavors of fermions with an approximate chiral SUl(N) x 
SUr(N) symmetry . The axial part of the above symmetry is broken spontaneously, when 
the fermion condensate (ipip) attains a nonzero vev. The low energy dynamics is then given 
in terms of the matrix field U(x) = exp(2i<f)(x) / /) in terms of the (A^ 2 — 1) Goldstone bosons 
4>(x). Here <ft(x) = <ft a (x)T a , a = 1,N 2 — 1 and T a 's are the generators of SU(N). Then 
the most general chirally invariant effective Lagrangian density with minimal number of 
derivatives given as [29] 

% = ^Tr(^Wt) (A.l) 
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with / as the corresponding 'pion' decay constant. An explicit symmetry breaking gives rise 
to finite masses of the Goldstone 'pions' and is given by 

C sb = ^^Tr (MC/t + C/Mt) (A.2) 

where M is the fermion mass matrix which breaks the chiral symmetry explicitly. In the 
following we shall also assume the the masses of the fermions to be the same rrif for all the 
flavors i.e M — rrifl in which case the symmetry breaking term can be written explicitly in 
terms of the PNGB fields as 

sb = N — - Trcos (— — ) ( A - 3 ) 

where summation over all the flavors is understood for the condensate. 

An explicit evaluation of the PNGB potential flA.3j) up to quadratic terms in the fields 



gives us the Gellmann-Oakes-Renner (GOR) relation relating the PNGB mass m with the 
explicit symmetry breaking scale rrif, 

™ 2 / 2 = -^(#). (A.4) 

To calculate the the temperature dependence of the condensate, we use the Feynman - 
Hellmann theorem, according to which the value of the condensate (V''0)t at finite temper- 
ature is related to the derivative of the free energy density with respect to the symmetry 
breaking parameter rrif , 

(^) T =(^) + J-n(T). (A.5) 

Assuming that the thermodynamic potential is dominated by the Goldstone modes, the 
free energy difference Cl(T) = f2(T) — Q(T = 0) is given as 

n(T) = J rfk ln(l - exp-^) (A.6) 



where, E = -\/k 2 + m 2 is the single PNGB energy. Next, we may use the fact that J^y = 

from the GOR relation Eq. flA.4j) and, eliminate the condensate (ipip) in favor of m using the 
same equation to obtain from Eq. (IA.5j) and Eq. (IA.6j) . 



(^)t = W)o (l - 2(N2 N \ ) , (A.7) 

where, 

1 f dk 1 1 T 2 



In the last step we have evaluated the integral in the chiral limit m «T [30 . 

Next, we proceed to calculate the leading temperature dependence of the masses of the 
Pseudo Goldstone bosons. This can be calculated in the same line as has been worked out 
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in Ref. 31] for two flavor case using Hartree approximation. The key quantity here is to 
calculate the pion self energy in the medium. To lowest order, the interaction part is given 
by collecting quartic order fields in the expansion of Eq.l lA.ll) and Eq. ( 1A.2I) and takes the 
form 



C 



int 



P [ 3 1 

^Tr(9 M 0)(^0 3 ) - ^Tr(d^ 3 )(d^) 
o o 



(A.9) 
(A.10) 




(a) 



(b) 



FIG. 5: PNGB self energy 



With the above interaction term, the contribution to the self energy is given by the tadpole 

The Hartree approximation 



diagram and the sun set diagram shown in Fig. 5(a) and 5(b) 



corresponds to evaluating these diagrams with the full single particle Greens function, rather 
than the bare ones. Explicitly within Hartree approximation, the PNGB self energy is then 
given as [31] 

2in \2N 2 -3) + (u 2 -k 2 + rh 2 ) 2A 



n(cu,k,T) 



t 



(A.ll) 



m v - . v- - - / 3 ^ 

In the above, the first term in RHS originates from the term proportional to m 2 in the 
Lagrangian the interaction term in Eq. (lA.lOj) . The other terms originate from the derivative 
couplings in the Lagrangian Eq.( lA.10|) acting respectively on the external PNGB propagator 
k 2 ) and on the PNGB loop (fh 2 ). fh is the in medium PNGB mass. Here we have 



ca- 



used the trace relation for the SU(N) generators Tr(T a T c T c T b ) = (N 2 - l)/(AN)5 ab and 
Tr(T a T c T b T c ) = —l/(4N)5 ab . further, we have here included the temperature dependent 
part of the propagator for the PNGB field and have discarded the vacuum co ntribution. The 
term t is the same as in Eq. (1A.8|) except the change E — » Cj = a/ {jh 2 + k 2 ), denoting the in 
medium PNGB dispersion relation with the effective mass fh. The effective mass satisfies 
the equation 



leading to 



m 2 



m 2 



m l + n(m,k = 0,T) 
3)ti 



2_ 

3N 

1 - 



4Ar 



-777. 



(A.12) 



(A.13) 
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further noting the fact that the GOR relation is also valid at finite temperature [3 lj, we obtain 
the temperature dependence of the the decay constant, using Eq. (IA.7l) and Eq.( lA.13]) .as 



/(T) = /(l-A^) (A.14) 



to the leading order in T/f. 
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